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INTRODUCTION 

In the past decades, many new discoveries in astron- 
omy that might have something to do with gravity are 
coming in. For example, stellar objects at the spiral arms 
of galaxies are rotating at faster speeds than that can be 
explained by the Keplerian motions. To overcome this 
difficulty, people assume that some extra matter, not vis- 
ible to us, is giving an extra pull on these stellar objects. 
Extra light deflections, as observed in the intergalactic 
gravitational lensing, are also ascribed to the existence of 
this extra matter. This is as known as the Dark Matter 
problem. An equally well known fact is that the Uni- 
verse is accelerating in its expansion. This is in contrast 
to our expectation that the Universe should be deceler- 
ating, unless some extra energy is kept pumping into the 
Universe. This is the so called the Dark Energy problem 

The expected Keplerian motions of stellar objects in 
the galactic spiral arms and the expected decelerating 
expansion of the Universe are something that arc pre- 
dicted by solving the Einstein Field Equation of General 
Relativity under the visibly observed mass-energy distri- 
butions in the spiral galaxies and in the Universe. The 
deviations from these expectations are attributed to, by 
most people, the existence of dark matter and dark en- 
ergy. But it could well happen that the General Theory 
of Relativity is not sophisticated enough to explain all 
these observed phenomena in gravity. Here we are tak- 
ing this point of view, and are proceeding to see if we can 
obtain an explanation of these anomalous astronomical 
phenomena by going beyond Einstein's General Theory 
of Relativity. 

A YANG-MILLS TYPE GAUGE THEORY OF 
GRAVITY 

There, by now, exist many modifications to Einstein's 
General Theory of Relativity so as to tackle gravity. 
Given the great success of the Yang-Mills way of de- 
scribing the electroweak and the strong interactions, we 
believe that the Yang-Mills way could also be a good 



approach to describe the gravitational interaction. 

The Yang-Mills type of action that describes gravita- 
tion will be taken to be the form of 

(1) 

where g^'^ is the metric and F^^, are the connections, and 
the Riemann curvature tensor is constructed out of the 
connections by 

This form of the gravitational action has presented it- 
self many times in the history of the development of the 
theory of gravity, but is, in fact, carrying very different 
information at each one of the presentations. The point 
of focus is on the relation between the metric and the 
connections. 

In its very early version, as proposed by Hermann 
Weyl [3, the connections that appear in the theory are 
nothing but the Christoffel symbols which are of first 
derivatives in the metric. This will result into a the- 
ory in which the metric is the only dynamical variable, 
and the variation with respect to the metric will give an 
equation of motion of higher order derivatives. It is well 
known that such a theory will possess runaway solutions. 

Later Yang also regarded the connections as the 
Christoffel symbols at the start, but varied the connec- 
tions instead in order to get the equation of motion. The 
final result is, again, an equation of higher derivatives in 
tlie metric. 

Stephenson put the anti-symmetric parts of the 
connections equal to zero, and regarded the symmetric 
parts of the connections and the metric as independent 
variables. And he obtained two equations of motion by 
varying both the metric and symmetric parts of the con- 
nections independently. 

On the other hand, some people identify the symmetric 
parts of the connections as the Christoffel symbols and 
regard the anti-symmetric parts and the metric as inde- 
pendent variables. Those people working on the Poincare 
Gauge Theory of Gravity are taking this point of view Q . 
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For us, we shall regard the full connections, both the 
symmetric and the anti-symmetric parts, as well as the 
metric as independent variables. The reason for us to 
take such a position is because the action in Eq. [1] with 
full symmetric and full anti-symmetric parts, is the Yang- 
Mills action for the GL{4: R) gauge group in the presence 
of a background metric . In this chosen theory of grav- 
ity, the connections are just the transformed GL{4 R) 
Yang-Mills vector potentials A^^, and can be taken as 
being independent of g^,^. 

The different choices of the content coded in the Ric- 
mann curvature tensor give different stories for physics. 
For example, for the Weyl theory, the metric is the only 
dynamical variable, and hence the action will contain ki- 
netic terms that have derivatives that are of orders higher 
than two. And when we look for the possible propagation 
modes in the theory, which can be obtained by looking at 
the inverse of the kinetic term, we will find that there will 
be propagators having the wrong signs, which will corre- 
spond to unphysical states called the ghosts or tachyons, 
and will end up into an unstable theory with the so called 
Ostrogradski instability . 

For our choice, the metric is a non-dynamical back- 
ground field 0, and the only dynamical variables arc 
the connections which obey an equation second order in 
space and time derivatives. The propagating modes are 
the 16 vector bosons and nothing else. And hence our 
chosen theory will contain no ghost and no Ostrogradski 
instability. 

THE TWO LEGITIMATE SOLUTIONS OF THE 
YANG-MILLS TYPE GAUGE THEORY OF 
GRAVITY 

The variation of the action given in Eq. [1] with respect 
to gsr will give the so called Stephenson Equation [3[ 

Her = R\gpR\^ - -^gerR^^^R^X^p = Tgr, (3) 

which as we can see, has an algebraic expression for the 
various components of the Riemann curvature tensor on 
the LHS of the equation. 

The variation of the same action with respect to the 
connections willgive the so called Stephenson-Kilmister- 
Yang equation [1, S, S] 

y p{T){^gRfP>^) = ^gSf\ (4) 

The Tqt and Sj^^ are respectively the metric energy- 
momentum tensor and the gauge current tensor of the 
source and the test object, coming from varying the mat- 
ter part of the action with respect to the metric and the 
connections. Vp(r) here denotes covariant differentiation 
with the connections F. 

If we choose to describe gravity with the above action, 
then it is mandatory for us to seek for all the possible 



simultaneous solutions to Eq. [3] and Eq. |4l These are 
complex equations, but we arc fortunate enough to lo- 
cate two spherically symmetric vacuum (where Tg-r ~ 
and S^^^ = 0) solutions in the literatures. They are 
both asymptotically Minkowskian and singular at r = 
(where the source of the gravity is supposed to sit at). In 
these two solutions, the anti-symmetric parts of the con- 
nections are zero (torsion free) and the symmetric parts 
happen to be the same as the Christoffel symbols calcu- 
lated from their respective metrics. These metrics are 
the Schwarzschild metric 

ds^ = (1 )dt^ - (1 y^dr^ - r^dn^, (5) 

r r 

and the Thompson-Pirani-Pavelle metric 

ds^ = {^l+':l^)-^dt^-{l + ^^r^dr^-r^dn\ (6) 
r r 

where GM and G'M' arc the integration constants for 
the solutions. 

In these two solutions, the metrics happen to be com- 
patible with their respective connections. But we have 
to bear in our mind that we have not assumed metric 
compatibility, as a priori, in our formulation of the the- 
ory. The existence of some solutions which are metric 
compatible does not mean that the theory is a higher 
derivative theory. 

The Thompson-Pirani-PavcUe solution was first dis- 
covered as a solution to the vacuum Eq. |4l with symmet- 
ric connections and with no reference to Eq. [3l It was 
later shown by Baekler, Yasskin, Ni, and Fairchild 
and by Hsu and Yeung [l3| that this solution satis- 
fies vacuum Eq. 2] for the full connections. It is also 
easy to see that these metrics will also satisfy the vac- 
uum Eq. [3] [^, [l3|. In fact, it was shown [3, that 
the Schwarzschild and the Thompson-Pirani-Pavelle so- 
lutions are the only two possible simultaneous solutions 
to vacuum Eq. [3] and vacuum Eq. |4] for spherical sym- 
metric situation under the compatibility ansatz. 

We don't know whether there exist solutions to our 
theory that have anti-symmetric components in the con- 
nections, or whether there exist solutions in which the 
metric is not compatible with the connections. But we 
shall assume that these yet undiscovered solutions will 
play no role in the discussions of the following physical 
phenomena. 

OUR INTERPRETATION OF THESE TWO 
LEGITIMATE SOLUTIONS 

Though the Schwarzschild solution has already found 
a lot of applications in the studies of various gravita- 
tional phenomena, the Thompson-Pirani-Pavelle metric 
was dismissed by its discoverers soon after its discov- 
ery because it failed to reproduce the classical tests that 
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are so successfully predicted by the Schwarzschild met- 
ric. And this leads, subsequently by many people, to the 
conclusion that the Yang-Mills type gauge theory is not 
a viable physical theory. 

Here we want to show that a suitable interpretation of 
the above two metrics can save this situation, and will 
reproduce nicely both the galactic rotation curves and 
the amounts of intergalactic lensing, if Nature follows 
our following postulate. 

We postulate that if Nature is going to make use of 
the Yang-Mills type theory of gravity, matter will then 
be endowed with either one of these two metrics. Mat- 
ter endowed with the first metric (denoted by g) will be 
called the regular matter, and matter endowed with the 
second metric (denoted by g') will be called the primed 
matter. Both the regular matter and the primed mat- 
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ter were produced during the creation of the Universe, 
though in different amounts, probably because of the dif- 
ference in the requirements of energy in producing them. 

There might possibly be other particles endowed with 
some yet undiscovered metric solutions. These parti- 
cles are assumed to be too heavy to be produced in an 
amount large enough to be noticeable in the present day 
astronomical observations. In other words, we will suffice 
ourselves by considering the physics played by these two 
gravitational copies of matter. 

Note that the energy- momentum tensor Tgrig) is a 
functional of the metric. When it is evaluated at it will 
be the energy-momentum tensor for the regular matter. 
And when it is evaluated at g' , it will be the energy- 
momentum tensor for the primed matter. So we have 



energy-momentum tensor of the regular matter = Tgrig), 
energy-momentum tensor of the primed matter = Tgrig'). 

I 



(7) 



It was explicitly pointed out by Stephenson Q that 
when we use the full connections (these connections are 
having both the symmetric and the anti-symmetric parts) 
as dynamical variables, then the Bianchi identities and 
vacuum Eq. |4] together will imply a local covariant con- 
servation law for the metric energy-momentum tensor 



V'iT)Hgr = v^r)re,(.9) = 0. 



(8) 



This is valid for any generic metric and connections which 
satisfy vacuum Eq. 21 In particular when the metric is g 
(corresponding connection T) or g' (corresponding con- 
nection r'), we will have two separate conservations laws 



v^f)r,.(5) : 

V'iT')Tgrig') 



0, 
= 0. 



(9) 



Eq.[7]and Eq.|9] together will mean that the regular mat- 
ter and the primed matter are separately conserved. 

Therefore, when we speak of a regular particle at x, we 
will mean an energy-momentum tensor that is localized 
at X and satisfies the first conservation law. And when we 
speak of a primed particle at x, we will mean an energy- 
momentum tensor localized at x and satisfies the second 
conservation law. 

Because a local covariant conservation law can be 
translated into an equation of motion for a test object 
which shares the same Tgr with the source but not af- 
fecting the metric , Eq. [7] and Eq. [9] will mean that 
the regular matter will move under the influence of the 
regular metric (which is generated by regular matter) 
while the primed matter will move under the primed met- 
ric (which is generated by the primed matter). In other 



words the regular matter will interact only with the regu- 
lar matter, and the primed matter will interact only with 
the primed matter, gravitationally. 

Other than gravitational interactions, the regular mat- 
ter and the primed matter are assumed to be identically 
the same in all other interactions. The respective ac- 
celeration produced by these metrics, when the speed of 
motion is small compared with the speed of light, can be 
calculated from 



2^ 



-dgt 



'00 



dr 



(10) 



and are respectively and -■^^;|^(1 + ^^^)^^- 

The G' and M' are introduced in order to make a 
parallel comparison between the Newtonian gravitational 
force and the new gravitational force. We will call G' the 
primed gravitational constant, and M' the primed gravi- 
tational mass. This new Gravitational force is attractive 
when G'M' is positive. 

There are some subtleties that we have to address care- 
fully in dealing with a geometric theory which has two 
gravitational copies of matter. Firstly, because the Grav- 
itational field equations are highly nonlinear, a rigorous 
solution for a source with both the regular matter M and 
the primed matter M' together is hard to find. However, 
we can expect that such a solution does exist. For such a 
solution, the metric will be reduced to the Schwazschild 
metric or the Thompson-Pirani-Pavelle metric when M' 
= or M = 0. Secondly, for a test object which itself is a 
mixture of the regular matter m and of the primed mat- 
ter m', the way how it reacts to the gravitational pull is 
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of high importance. According to Eq. [7] and Eq. [9l each 
component will move in accordance with their own corre- 
sponding gravitational pulls. The combined gravitational 
force will be 



F 



GM 



G'M' 



(1 



G'M' , 



(11) 



and then we will arrive at a relationship holding between 
the rotation speed v with the distance r, 



GM 



G'M' 



r + G'M' m + m'' 



(14) 



The above force equation can also be understood by 
noting that when the test object is far from the source, 
it is in the weak field limit where the connections will 
be a superposition of the connections from M and from 
M', and the local covariant conservation will become, 
symbolically 



dT + (T + r')T = 0. 



(12) 



When integrating the above equation over a small space- 
time region surrounding the regular matter in the test 
object, it will give the regular gravitational force. And 
when integrating the above equation over a small region 
surrounding the primed matter, it will give the new grav- 
itational force. Note that even though the regular matter 
and the primed matter may be very close together, they 
nevertheless occupy different spacetime points and hence 
their energy-momenta will not overlap, and hence the 
above reasoning work well. 

When the regular matter and the primed matter are 
bound together by non-gravitational force, and move to- 
gether in circular orbit, the centrifugal force will balance 
the combined gravitational forces 



(m + m )- 



GM 



G'M' 



(1 



G'M', 



THE YANG-MILLS TYPE GAUGE THEORY OF 

GRAVITY PREDICTS THE UNIVERSAL 
ROTATION CURVES FOR SPIRAL GALAXIES 



There is an immediate application of Eq.[T4]to describe 
the kinematics of a stellar object which is made up the 
regular matter m and the primed matter m' , and is mov- 
ing at a distance r from the center in a spiral galaxy. 
This spiral galaxy whose disk is assumed to be composed 
of solely of regular matter of mass M, and whose halo is 
assumed to be composed of a uniformly distributed mat- 
ter with a contamination of the primed matter of mass 
density p' . Because 



P 



(15) 



(13) 



M'^^r^- 
3 

where M' is the total primed mass inside r, and because 
the arm structure of the regular matter presents its con- 
tribution in modified Bessel functions form [l^l , the total 
contributions to by the spiral arms and the halo will 
sum up to 



m + nv r + G'-^p'r-^ m + m' 

47rGI]oV[/o(y)A'o(y) - Ii{y)K,{y)] x (— ^) + ^tn^^—^^^ (^6) 

m + m' r + G*r^ m + m' 



3 



m 



where /i, Ki arc modified Bessel functions of the first 
and second kinds. The is the speed coming from the 
Newtonian force of the disk with y = ^ [13] ■ 

We will see immediately that this agrees exactly with 
the empirical formula given by Salucci et al [13| who 
analyzed a large number of spiral galaxies and drew up 
a formula to describe their rotation curves. We have 
also extracted the values of G* (which is the product of 
the primed gravitational constant and the primed mat- 
ter density in the halo) and m'/m (which is the ratio 
of the primed mass density to the regular mass density 
in the halo) by fitting some well-known galaxies [IJ- 
18j . and it is amazing to find that a universal value of 



G* w 10 ^ kpc ^ and a universal ratio of — « 2x10 ^ 
fit very well with the observed results when r ranges from 
3 kpc to 30 kpc. Note also that the value of GEq and h 
are more or less the same as those observed. The results 
are shown in Fig. [T] and Table HI 

We also have to check for the influence of this diffuse 
halo medium on the motion of the planets in our solar 
system. With the values of the G* and m'/m given in 
the above, the observed rotation speeds of the planets fit 
well with the predictions of Eq. [161 Figure [2] gives the 
planetary motions for various values of G* and m'/m. 

The alert reader may find that, in the above discus- 
sions, we have already made the assumption that the 
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FIG. 1: The predicted relationship between the galactic rotation speed v and the distance r from a combined influence of the 
Newtonian force and the new gravitational force: [(a)]rhe Milky Way, [(bjjMGC 3198, [(c)]N'GC 2403 andp)]MGC 6503. 





The Milky Way 


NGC 3198 


NGC 2403 


NGC 6503 


GSo [km''s-''kpc-'] 


6.8 X 10^ 


2.8 X 10^ 


2.1 X 10^ 


2.8 X 10^ 


h [kpc] 


2.0 


2.63 


2.05 


1.72 


m. 


5.0 X 10-2 


9.2 X 10-3 


1.4 X 10-2 


1.3 X 10-2 


2.3 X 10^9 


2.2 X 10-9 


2.0 X 10-9 


1.4 X 10-9 



TABLE I: The fitting parameters for some galaxies. 



stellar objects in the galaxies are composed solely of reg- 
ular matter. This assumption can be understood in the 
following way. The primed matter always respond to the 
primed gravitational pull with a high rotation speed even 
when they are bound with some regular matter and hence 
are far harder for them to condense gravitationally into 
a star. 



Though the stars contain no primed matter, the stars 
that are rotating at the outskirt of a spiral galaxy are, in 
fact, embedded in pockets of halo media and are rotating 
around the galaxy all together. 

The amount of regular matter in the halo is small when 
compared with that in the galactic bulge and spiral arms, 
and is thus neglected in the above discussions. 
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FIG. 2: The orbital speeds of the planets in the solar system 
predicted by Eq. [TS]with different m'/m values. 2 x 10~', 
1 X 10~* and 1 x 10^^ are for red, blue and green points 
respectively. 



THE YANG-MILLS TYPE GAUGE THEORY OF 
GRAVITY GIVES THE RIGHT AMOUNTS OF 
INTERGALACTIC GRAVITATIONAL LENSING 

Next, let us turn to see what the primed matter does 
in explaining the large light deflections that arc observed 
in intergalactic gravitational lensing. We shall take the 
galaxy cluster Abell 1689 as our illustration. We shall re- 
gard Abell 1689 as a cluster consisting of galaxies which 
are carrying their own individual halos with them. And 
this collection of galactic halos forms the halo of the clus- 
ter. Since the galactic halos are always regarded as hav- 
ing a size of the order of 30 kpc, the size of the cluster 
will be very close to its halo size which is taken to be 300 
kpc. 

The Azimuthal angle swept by the light, when it travels 
from point R to the point of closest approach rg, under 
the influence of gravity described by the metric 

ds2 ^ B{r)dt^ - A{r)dr^ - r^dil^, (17) 



is given by [l9| 

Jro ^0 B(r) r 

(18) 

In the case of a Thompson-Pirani-Pavelle metric for a 
point source of primed mass A/', the angle swept is 

Aip 



(rp + G'M')dr 

(r + G'M')[(r + G'M'f - (ro + G'M'Y]^' 

(19) 

which can be readily integrated to give 



sec 



(20) 



with P = ^TXp. 

We will immediately notice that the angle change will 
be 7r/2 when R goes to inflnity. That means no deflec- 
tion for the light by a point source of the primed matter 
when it comes from infinity and then goes back to infinity. 
Things will be different when we deal with a distributed 
source of the primed matter as we are going to show in 
the following. 

Let R be the radius of the halo of Abell 1689 and tq be 
the closest approach from the cluster center. The light 
will see a point source of the primed matter of constant 
mass M' when it is moving beyond the halo. And it will 
see a point source of diminishing mass (i.e. increasing /?) 
when it enters the halo because it sees only the mass that 
lies inside r. 

We claim that the deflection by a point source of mass 
M' is smaller than that by a uniformly distributed halo 
which has a total mass of M', if the light penetrates into 
the halo during some time in its journey. The above claim 
is obvious by noting that 



d{Aip) 
dp 



> 0. 



(21) 



So if the light penetrates into the cluster halo, the total 
Azimuthal angle change will be greater than 7r/2, and we 
will see light deflectng towards the center of the cluster. 

The actual angle swept, when light comes from the 
infinity, enters the halo at R and reaches the point of 
closest approach at rg is given by 



+ (A(y9|halo — A</j|point source)- 



(22) 



The second term in the bracket in Eq. [21] comes from 
the fact that we have to take out a 7r/2 in calculating the 
angle of deflection. 

The quantity inside the bracket of Eq. [22] can be re- 
garded as a change S{Alp) in Aip due to a change SM' in 
Af, and we get 



Aiy9|halo - 

= S{A^) 
_ djAip) 
dM' 
d , 



A^\ 



point source 



(23) 



dM 
3 

71 



1 



{R-rQ)iG*~^R 



]SM' 



A plug in the data of i? = 300 kpc, (i? — Tq) = 30 kpc 
and G* = 10"^ kpc"^, will give us an angle of deflection, 
which is 2d{Aip), a value of 4 x 10^'^, which is many times 
of that expected from General Relativity. 

Note that the value of G* that we used in calculat- 
ing the Abell 1689 light deflection comes from the curve 
flttings of galactic rotations. Again there seems to be a 
universal value for G*. as it should be. 
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We should also note that the sun and the planets in the 
solar system aren't contaminated with the primed matter 
as we have explained in the above. The fundamental tests 
on General Relativity will hence remain intact. 



THE PRESENCE OF A PRIMORDIAL TORSION 
EXPLAINS THE ACCELERATING EXPANSION 
OF THE UNIVERSE 

There is another nice feature of the Yang-Mills type 
gauge theory of gravity, when we use it to study the Uni- 
verse as a whole. This theory was shown to admit a 
cosmological solution of the form 2^ 2l| 



Hehl and Professor Jim Nester for valuable discussions. 



with primordial torsion compoments, 

Tioi = ^202 ~ Tsoa = ^. 



(24) 



(25) 



The po and ^, with ^ > are integration constants 
arising from the integration of the equation of motion. 

We can interpret this solution as representing an ex- 
panding and accelerating Universe when the influence of 
gravity dominates over the influence of matter and radi- 
ation. 

The role played by the primordial torsion is crucial 
here: the stretching on the Universe by the metric is 
compensated by the twisting by the torsion. And from 
the metrical point of view, we look like living in a Uni- 
verse with a cosmological constant ^. 

And interesting enough, our torsion selects the spa- 
tially flat metric {k — 0) as the only accompanying met- 
ric [2^, [2l| . The spatially flat geometry of the Universe 
is confirmed by WMAP. 



CONCLUSION 

The Yang-Mills type gauge theory of gravity has a 
richer structure than that of Einsteins General Theory 
of Relativity. Its richer number of solutions, both in the 
absence and in the presence of torsion allow us to describe 
more physical phenomena with it. The recently observed 
astronomical phenomena seems to show that Nature is 
enjoying the full use of the Yang-Mills type gauge theory 
of gravity. 
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